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SUMMARY 


Exact  and  approximate  methods  based  upon  linearized  supersonic  flow 
theory  have  been  developed  for  the  calculation  of  the  velocity  potential 
and  the  downwash  from  thin  wings  of  arbitrary  plan  form.  Particular 
attention  is  given  to  the  evaluation  of  the  downwash  in  the  plane  of  the 
wing.  The  applicability  of  the  method  inherently  depends  upon  a  knowledge 
of  the  load  distribution  over  the  plan  form  of  the  wing.  General  expres¬ 
sions  for  the  velocity  potential  and  downwash  have  been  derived.  Simple 
modifications  of  these  expressions  produced  formulas  for  the  velocity 
potential  and  downwash  from  arbitrary  curved  lifting  lines. 

A  complete  development  of  all  formulas,  starting  with  the  basic 
solution  of  the  linearized  potential  equation  for  supersonic  flow,  is 
given.  Although  the  paper  contains  many  new  results,  some  of  the  results 
presented  have  been  obtained  by  other  methods  and  are  given  here  solely 
for  completeness.  The  general  formulation  of  the  downwash  equations  can 
easily  be  used  in  finding  exact  and  approximate  expressions  for  the  other 
velocity  components. 

The  results  of  the  theoretical  development  are  used  to  determine  the  * 
downwash  from  a  pitching  rectangular  wing  and  to  determine  the  expression 
for  the  local  angle  of  attack  necessary  to  give  a  specified  load  distri¬ 
bution.  Comparisons  of  the  exact  and  approximate  values  of  the  downwash 
for  several  lifting-line  configurations  are  also  presented. 


INTRODUCTION 


The  calculation  of  the  downwash  fields  induced  by  thin  wings  at 
supersonic  speeds  is  necessary  in  order  to  evaluate  accurately  the  aero¬ 
dynamic  load  distribution  over  the  tail,  an  important  consideration  for 
structural  and  stability  calculations.  Available  methods,  including  the 
present,  are  based  on  the  linearized -time -independent  flow  considerations. 
These  methods,  in  general,  utilize  the  following  well-known  concepts: 
conical  flows,  potential  doublets,  vortices,  and  pressure  doublets.  The 
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conical-flow  method  is  used  in  reference  1  to  calculate  the  downwash 
flow  fields  in  the  plane  of  the  wing  for  several  plan  forms  including 
the  rectangular  wing.  The  potential -doublet  method  is  presented  in  ref¬ 
erence  2  and  is  used  to  find  the  downwash  behind  a  triangular  wing  over 
a  range  of  Mach  numbers  in  reference  3*  The  vortex  method  is  developed 
for  supersonic -flow  applications  in  references  4  and  5  &nd  is  used  in 
reference  6  to  find  the  downwash  due  to  thin  wings  approximated  by  lifting 
lines.  The  method  used  in  reference  7  was  the  potential-doublet  method; 
however,  by  integration  by  parts  Ward  obtained  an  expression  for  the 
velocity  potential  in  space  which  agrees  with  the  expression  for  the 
velocity  potential  in  space  determined  herein. 

The  present  method  is  essentially  a  development  of  the  pressure- 
doublet  method  for  the  calculation  of  downwash  flow  fields.  This  method 
leads  easily  to  expressions  for  the  velocity  potential  and  the  downwash 
from  arbitrary  curved  lifting  lines.  These  expressions  in  turn  lead  to 
approximate  expressions  for  the  downwash  from  lifting  lines.  These 
approximate  expressions  have  many  computational  advantages. 

The  method  described  herein  has  certain  advantages  in  that  most  of 
the  important  results  obtained  by  other  methods  can  be  obtained  without 
difficulty  by  using  the  present  approach.  Integrations  are  performed 
only  on  the  plan  form;  whereas  other  methods,  excluding  the  conical-flow 
method,  generally  lead  to  integrations  on  the  plan  form  and  over  the 
wake.  An  attempt  was  made  to  present  a  fairly  complete  development  of 
the  pres  sure -doublet  method  and  at  the  same  time  present  results  which 
have  not  been  obtained  by  other  methods. 

The  accuracy  of  the  developed  approximate  formulas  is  brought  out 
in  the  evaluation  of  a  rectangular  pitching  wing  by  use  of  exact  (lin¬ 
earized)  and  approximate  expressions.  Comparisons  are  made  between 
exact  calculations  of  the  downwash  from  certain  lifting  lines  and 
approximate  expressions.  Other  applications  include  some  simple  deriv¬ 
ations  of  known  results  plus  an  expression  for  the  local  angle  of  attack 
necessary  to  give  a  specified  load  distribution. 


SYMBOLS 


chord 
semi span 

X]_  limits  of  integration 
pressure  coefficient  ^Pressur^—pV^ 
denotes  finite  part  of  integral 


c 

b/2 

l1*l2 

P 

f 
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f(yl) 

equation  of  an  arbitrary  line 

i 

variable  index 

k 

a  constant 

hi,h2 

Yl  limits  of  integration 

M 

free-stream  Mach  number 

m 

slope  of  line 

n.n-Ljng 

limits  of  summation 

q 

*  . . 

angular  velocity  about  y-axis 

R-  =  yx2  -  P2(y2  +  Z2) 

R  =  /x2  -  P2(y2  +  z2) 

R(X)  =\J(\  -  x-jJ2  -  P2(y 

2*z2) 

V 

free-stream  velocity 

X  =  x  -  xlf  Y  =  y  -  yi 

>  z  “  z  “  zl>  xi  =  x  -  xi*  Yi  =  y-yi 

x»y>z>x1,y1,z1 

Cartesian  coordinates 

5xl>&yi 

increments  in  and  y^,  respectively 

a 

angle  of  attack 

p  =  \/m2  -  1 

r 

circulation 

e 

a  small  positive  number 

X 

an  auxiliary  variable 

P 

free-stream  density 

a 

local  angle  of  attack  :  *  ' 
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T 


T* 


* 

S0 


area  of  plan  form  in  forward  Mach  cone  from 
point  (x,y,z) 

area  of  plan  form  and  wake  in  forward  Mach  cone 
from  point  (x,y,z) 

perturbation  velocity  potential 

increment  in  velocity  potential  at  point 
(x,y,z)  due  to  an  elementary  lifting  area 
at  point  (x-^y-^o) 

difference  in  partial  derivative  of  0  with 
respect  to  xi  ^X]_u  "  $x±l) 

value  of  potential  difference  at  trailing  edge 
(0-PEu  "  ^TEl) 

vertical  perturbation  velocity  at  point  (x^y^z) 
due  to  regions  I,  II,  and  III,  respectively 


line  integral  around  area  of  plan  form  in 
forward  Mach  cone  from  point  (x,y,0) 


* 


Subscripts : 
u,  Z 
x,y,z 

TE 


upper  and  lower  surfaces,  respectively 

partial  derivatives  with  respect  to  x,  y, 
and  z,  respectively 

trailing  edge 


THEORY 


The  theoretical  development  is  divided  into  three  parts:  First, 
general  formulas  are  derived  for  the  potential  and  the  upwash  in  space 
due  to  a  thin  lifting  surface;  the  pressure  distribution  on  the  surface 
is  assumed  to  be  known.  Second,  because  the  formulas  in  the  first  part 
are  in  rn«my  cases  difficult  to  evaluate,  the  expressions  for  the  thin 
wing  are  used  to  develop  formulas  for  the  potential  or  upwash  m  space 
due  to  lifting  lines.  Third,  because  in  some  cases  the  expressions  for 
lifting  lines  are  difficult  to  evaluate,  approximate  expressions  are  also 
derived  for  the  upwash  from  lifting  lines. 
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Lifting  Surfaces 

Potential  In  space.-  The  partial-differential  equation  satisfied  by 
the  perturbation  velocity  potential  in  supersonic  flow  is 


p2<fg  _  _  &£  =  0  (!) 

8z^ 

A  general  solution  of  equation  (l)  is  given  in  reference  8.  (The  boundary 
conditions  for  airfoils  are  given  in  reference  9*)  For  the  purposes  of  the 
present  paper  this  solution  may  be  written  in  the  form 


^ w)  -  2i  jT't,  <(<u  '  dyi  (2) 

where  0U  and  0^  are  'tiie  values  of  0  on  the  upper  and  lower  sides  of 
the  surface  z  =  z The  finite  part  of  the  integral  in  equation  (2) 
must  be  taken  as  indicated  by  the  symbol  f .  The  double  integrals  that 
arise  will  be  dealt  with  according  to  the  methods  for  finding  the  finite 
part  of  the  multiple  integrals  given  by  Hadamard  in  reference  8.  The 
area  of  integration  Tf  is  the  area  of  the  z  =  z^  plane  for  which 

Xl  <  x  -  3  V^2~Tz2 

Equation  (2)  remains  valid  when  the  velocity  potential  is  replaced 
by  its  partial  derivative  with  respect  to  the  free-stream  direction. 

Since  the  pressure  is  directly  proportional  to  0X,  replacing  0  by  0X 
in  equation  (2)  has  the  effect  of  introducing  a  potential  which  is 
directly  proportional  to  the  acceleration  potential. 

Since 


0(x, y,z) 


0x(x,y,z)d\ 


it  follows  from  equation  (2)  that 


(3) 


0(x,y,z). 


(<v  -  dyi  ax 


(M 
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vhere 


R'(X)  =y(X  - >i)2  -  32(y2  +  Z2) 


The  difference  in  0Xl,  A0Xl>  is  zero  except  on  the  plan  formj  therefore, 
the  region  of  integration  T  is  over  the  plan  form  only. 

It  is  proved  in  reference  k  that  the  order  of  integration  of  the 
finite  part  of  a  multiple  integral  can  be  changed.  Equation  (4)  for 
zi  =  0  can  therefore  be  expressed  in  the  form 


0(x,y,z)  =  -L  JJ  (0X]LU  -  0xl2)  J  ^  dX  dxi  dyx 


The  integral 


f  0x 

r 


must  he  integrated  from  the  aftercone  emanating  from  point  (xi,yi,0) 
to  the  point  (x,y,z)  because  a  disturbance  at  point  (x^,y^,0)  only 
affects  points  in  the  aftercone  from  the  point  (xi^yi^O).  The  value  of 


the  preceding  integral  is 


(Y2  +  z2)r 


Equation  (5)  now  becomes 


0(x,y,z)  =  — 


2n  JJt  (Y2  +  Z2)R 


dx-L  dy± 


In  general,  equation  (6)  is  a  finite  integral;  therefore,  the  finite- 
parts  symbol  has  been  dropped.  The  expression  for  the  velocity  potential 
given  by  equation  (6)  was  obtained  in  reference  7  by  a  different  approach 
as  indicated  in  the  introduction.  Equation  (6)  may  also  be  readily  obtained 
from  Volterra’s  solution  (reference  10). 
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Upvash  in  space.-  Perhaps  the  most  frequently  desired  quantity  is 
the  upwash  jflz.  It  is  shown  in  the  appendix  that  equation  (6)  may  be 
differentiated  with  respect  to  z  under  the  integral  signs  without  con¬ 
sideration  of  the  variable  limits.  The  result  of  differentiating  both 
sides  of  equation  (6)  with  respect  to  z  is 


z2)r2  +  02z2(y2  +  z2  jlx 
- - - 2L.  dx-, 

(y2  +  z2)%3 


dyi 


(7) 


An  expression  for  the  sidewash  due  to  a  discontinuous  pressure  sheet  can 
be  obtained  in  a  manner  similar  to  that  used  in  finding  equation  (7). 

Upwash  in  the  wake . -  In  the  z  =  0  plane  for  points  not  on  the  plan 
form  or  in  the  wake,  equation  (7)  reduces  to  the  form 


1.  rr  A0Xl  X 

2  it  Jo  f  Y2R 


dx- 


The  restrictions  on  equation  (8)  can  easily  be  understood  by  comparing 
it  with  equation  (6).  Equation  (6)  is  an  expression  for  the  potential 
in  space  due  to  a  discontinuous  pressure  sheet.  For  points  not  on  the 
plan  form  or  in  the  wake,  it  can  be  seen  from  equations  (6)  and  (8)  that 


0z(x>y>°)  =  iig  \  0(x,y,z) 


The  preceding  expression  shows  that  in  order  for  the  right  side  of  the 
expression  to  be  finite  the  potential  must  approach  zero  as  z  approaches 
zero.  It  is  well-known  that  the  potential  is  zero  in  the  z  =  0  plane 
for  zero-thickness  wings  except  on  the  plan  form  or  in  the  wake  where  the 
potential  has  a  finite  discontinuity  across  this  plane.  Because  0Z  is 
general  finite  and  since  the  potential  has  a  finite  discontinuity 
across  the  z  =  0  plane  in  the  wake  and  on  the  plan  form,  equation  (8) 
cannot  hold  on  the  plan  form  or  in  the  wake. 

Equation  (7)  may  be  used  to  find  0Z  for  points  on  the  plan  form 
or  in  the  wake  by  performing  the  integrations  for  an  arbitrary  z  and 
then  setting  z  equal  to  zero.  Generally  the  integrations  for  an  arbi¬ 
trary  z  are  very  difficult.  Both  integrations  of  equation  (7)  must  be 
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performed  for  points  on  the  plan  form;  whereas,  only  the  integration  with 
respect  to  yp  is  needed  for  points  in  the  wake.  The  reason  that  one  or 

both  integrations  must  he  performed  before  taking  the  limit  of  0Z  as  z 

approaches  zero  is  that  the  integrand  tends  to  infinity  along  the  line 
y-^  =  y  as  z  approaches  zero.  The  line  yp  =  y  must  therefore  be 
removed  from  the  area  of  integration.  This  line  may  be  removed  by 
removing  a  narrow  strip,  of  width  6 ,  on  each  side  of  the  line  yi  =  y 
from  the  area  of  integration.  (See  fig.  1(a).)  The  integrations  must  be 
performed  until  the  width  e  is  removed  from  the  limits  of  integration. 
The  limit  of  0Z  as  z  approaches  zero  may  then  be  taken.  The  limit 
of  0Z  as  z  approaches  zero  cannot  he  taken  until  the  limit  of  0Z 

as  6  approaches  zero  is  taken  unless  the  effect  of  the  strip  is  con¬ 

sidered.  Figure  1(b)  shows  the  area  of  integration  necessary  to  evaluate 
a  point  in  the  wake.  It  can  be  seen  from  this  figure  that  the  xp  limits 
are  independent  of  €;  therefore,  for  points  in  the  wake  the  xp  inte¬ 
gration  need  not  be  performed  before  the  limits  are  taken. 

The  calculation  of  0Z  for  points  on  the  plan  form  or  in  the  wake 
can  be  considerably  simplified  by  using  a  combination  of  equation  (7)  and 
equation  (8).  For  points  in  the  wake,  divide  the  area  of  integration  as 
shown  in  figure  2.  Each  region  of  the  integration  may  be  considered  as 
an  independent  wing  with  a  local  angle  of  attack  necessary  to  give  the 
load  distribution  of  the  region.  Since  regions  I  and  III  can  be  con¬ 
sidered  as  being  two  different  plan  forms,  equation  (8)  can  be  applied 
to  regions  I  and  III,  but  equation  (7)  must  be  applied  to  region  II. 

The  value  of  e  is  chosen  small;  therefore,  in  region  II  the  variation 
of  A0X  with  yp  can  be  neglected  in  the  first  approximation  provided 
the  pressure  is  continuous  over  the  plan  form.  Pressure  discontinuities 
may  be  dealt  with  by  subdividing  the  regions  of  integration. 

The  contribution  of  regions  I  and  III  to  the  vertical  perturbation 
velocity  when  z  =  0  is 


z'l+III  2  it 


A0xi  X 


dxx  dyp 


I+III 


For  a  point  above  the  wake  the  effect  of  region  II  is  approximately 


'z'll  ~  2 it 


L0  n y+e 


Lp  Jy-e 


A0xp(*i>y>Q)x[(Y2  -  z2)r2  +  p2z2(y 

^Y2  +  z2)2R3 


2  2 
^  +  zc 


dxx  dy1 


(10) 
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where  A0x(xi,y,O)  is  the  difference  in  0X  across  the  lifting  surface 
along  the  line  yj_  =  y.  The  xj_  limits  are  denoted  by  Lj_  and  L2 • 
After  integration  with  respect  to  equation  (10)  becomes 


<^>ii  *kjL 


L2  A0Xl(xi,y,O)XY(R2  -  p2z2) 
Ll  (x2  -  P2z2)(y2  +  z2)r 


y+e 


y-e 


dx- 


and  for  z  =  0 


<d  \  ~  1  rL2  A^xi(xi^^0)  - p2e^ 

zn~"^JLl  X 


dXf 


(id 


For  points  in  the  wake  the  resultant  vertical  perturbation  velocity 
is  then  given  by  the  sum  of  equations  (9)  and  (ll): 


0z(x>y>°) 


'XI 


i+iii  y2r 


dxx  dy! 


i  [h 

rte  Jl± 


l2  A(2(Xl(xi,y,0)  \jx2  -  p2e2 


dx- 


(12) 


Equation  (12)  is  approximate  when  €  is  not  zero;  however,  the  accuracy 
of  the  approximation  increases  as  6  approaches  zero  and  of  course  in 
the  limit  the  equation  is  exact.  The  exact  vertical  perturbation  velocity 
is  then  given  by 


0z(x,y,O)  =  g1^ 


JL 


A0xi  X 


2jr  JJ I+III  Y2R 


dyi  dxx  - 


1  Pl2  A0Xl(xi,y,O)\/X^  - 


—  r  ‘ 

Jte  JLl 


,2,2 


dx- 


(13) 
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Equation  (13)  can  be  expressed  in  the  form 


0z(x>y>°)  = 


”&JL 


’xiyi 


R 


XY 


dxp  dy-L 


•*£ 


!xi \I 


X2  -  p2y2 


XY 


dx-L  (14) 


where  the  second  integral  is  a  line  integral  around  the  region  of  the 
plan  form  affecting  the  point  (x,y,0).  Equation  (l4)  was  obtained  from 
equation  (13)  by  integrating  the  first  term  by  parts  with  respect  to  ylv 

A  consideration  of  the  previous  division  of  the  plan  form  into 
regions  with  regard  to  vortex  distributions  is  helpful  in  understanding 
the  physical  meaning  of  the  preceding  manipulations. 

Equation  (10)  is  the  expression  for  the  upwash  from  a  series  of 
horseshoe  vortices  distributed  over  region  II.  The  y^  component  of  the 

vortex  strength  is  proportional  to  A0x(xq,y,O) .  The  xq  component  of 

the  vortex  strength  is  zero  except  along  y  -  e  and  y  +  e.  (See  fig.  3*) 
The  spanwise  distribution  of  the  trailing  vortices  consists  of  only 
two  vortices  of  finite  strength  A0(x,y,O)  located  along 


Yl  =  y  -  e 


and  along 


yq  =  y  +  6 


■* 


Equation  (9)  can  be  considered  as  the  expression  for  the  vertical 
perturbation  velocity  due  to  a  system  of  vortices.  The  integral  over 
region  I  represents  the  effect  of  the  bound  and  trailing  vortices 
associated  with  region  I  plus  a  finite  vortex  (see  fig.  3)  along 


yq  =  y  -  e 

This  finite  vortex  is  the  sum  of  the  y-components  of  thd  vortex  strength 
along 


yq  =  y  -  e 
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Since  the  divergence  of  a  vortex  field  is  zero,  it  is  not  surprising 
that  a  vortex  of  finite  strength  exists  along  this  line.  The  integral 
over  region  III  represents  the  effects  of  the  hound  and  trailing  vortices 
associated  with  region  III  plus  a  vortex  of  finite  strength  along 


yi  =  y  +  € 


Equation  (9)  gives  the  effect  on  the  field  point  of  the  vortex  distri¬ 
bution  on  regions  I  and  III,  the  trailing  vortices  associated  with  the 
vortices  on  regions  I  and  III,  and  two  vortices  of  finite  strength 
located  a  distance  €  on  either  side  of  the  line  yi  =  y.  For  any  small 
value  of  e  the  strengths  of  the  two  finite  vortices  are  almost  equal 
and,  since  the  field  point  lies  midway  between  these  finite  vortices,  it 
follows  that  equation  (9)  should  tend  to  infinity  as  €  approaches  zero. 
The  infinity  arising  from  equation  (9)  as  €  0  is  exactly  canceled  by 

the  infinity  arising  from  equation  (ll)  as  e  •>  0.  A  finite  downwash 
then  results  at  the  point  (x,y,0)  within  the  wake  region  provided  that 
the  spanwise  derivative  of  the  load  distribution  is  continuous  at  y. 


Upwash  on  the  plan  form.-  In  order  to  simplify  the  calculation  of  0Z 
on  the  plan  form,  consider  an  arbitrary  point  slightly  above  the  plan 
form  and  divide  the  area  of  integration  into  regions  as  shown  in  figure  4. 
Region  IV  is  small;  therefore,  the  value  of  A0X-^  is  approximately 

constant  in  region  IV.  For  points  slightly  above  the  plan  form  the 
contribution  of  region  IV  to  the  vertical  perturbation  velocity  is 
approximately 


(0z) 


IV 


2« 


7 r  d 

i*2  -  *2) 

Ir2  +  (32z2| 

[y2  +  z2)] 

IJ  iv 

•(Y2  +  Z2)2r3 

dxx  dy-L  (15) 


The  result  of  performing  the  integrations  in  equation  (15)  is 


(0z)xv  ~  -f-  A0x(x,y,O)  (16) 

Actually  the  point  at  the  apex  of  the  hyperbola  is  a  singular  point  and 
must  be  removed  from  the  region  of  integration  by  a  method  such  as  is  given 
in  reference  8,  pages  1V7  to  150.  The  integrals  over  regions  I,  II,  and  III 
for  points  slightly  above  the  plan  form  are  given  approximately,  by  equa¬ 
tion  (l4) .  The  vertical  perturbation  velocity  fz  is  then  approximately 
given  by  (equation  (l4)  plus  equation  (l6)): 
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A0X]y]_ 


2  )t  JJT  XY 


dXx  dy-L 


.l£ 

2  it  Jm 


A0Xi  R 


MC  XY 


dX]_  -  |  A0x(x,y,O)  (17) 


Equation  (17)  is  independent  of  e ;  therefore,  it  is  exact  for  points  on 
the  plan  form.  The  result  given  by  equation  (17)  was  also  obtained 
independently  by  Dr.  A.  Busemann  of  the  Langley  Laboratory  in  an  unpub¬ 
lished  analysis  using  a  different  approach.  <■ 


Lifting  Lines 

Integrations  that  arise  in  evaluating  the  potential  or  the  downwash 
from  thin  wings  are  in  most  cases  very  difficult  to  perform.  For  most 
downwash  problems  reference  6  indicates  that  a  lifting  line  can  be  used 
as  a  very  good  approximation.  In  this  reference,  formulas  are  developed 
for  the  upwash  in  space  due  to  lifting  lines.  In  the  present  paper, 
formulas  are  developed  for  the  velocity  potential  and  the  upwash  in  space 
due  to  lifting  lines  by  using  the  same  approach  that  was  used  to  develop 
the  previous  formulas  for  the  thin  wings .  The  results  for  the  upwash 
agree  with  the  corresponding  results  of  reference  6. 

Potential  due  to  a  lifting  line.-  The  potential  or  upwash  due  to  an 
arbitrary  curved  lifting  line  can  be  found  from  the  preceding  results. . 
The  following  expression  for  the  infinitesimal  increment  in  the  potential 
due  to  an  elementary  lifting  area  can  be  obtained  from  equation  (6): 


80(x,y,z)  =  A0xi  zX 
6xl  8yl  2it(Y2  +  z2)R 

The  potential  from  an  element  of  a  lifting  line  is  then 


(18) 


d0(x,y,z) 


A0xn  5xl  zX  5Yl 
lim  - - 

Sx^y^O  2it(Y2  +  z2)R 


(19) 


where  the  product  A0Xl  Sxp  is  held  constant  as  5x-j_  approaches  zero. 

The  product  A0Xl  is  the  difference  in  potential  because  A0Xl  is 

constant  in  the  x-direction.  Since  A0  '  is  also  the  circulation,  equa¬ 
tion  (19)  can  be  written 
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d0(x,y,z)  =  - - dy-L  (20) 

2Jt(Y2  +  z2)R 

where  V  denotes  the  circulation  at  point  (x]_,yq). 

Let  the  equation  of  the  curve  denoting  an  arbitrary  lifting  line  be 


*i  =  f(yi) 


(21) 


Substituting  this  function  into  equation  (20)  and  integrating  yields  the 
equation 


0(x,y,z) 


i_  rh2  z(x  -  f)r(yi)  d 

2lt  (Y2  +  z 2)\/(x  -  f)2  -  (32(Y2  +  z2)  yi 


(22) 


Upwash  due  to  a  lifting  line . -  The  derivative  of  equation  (22)  with 
respect  to  z  is 


2% 


T*l 

[(x-  f)2-  (32Y2 

Y2- 

z2 

[(x  -  f)2-  p2(Y2  +  2z2)] 

}(x-f)r 

■4 

(Y2+  : 

z2)2[(x 

-  f  )2 -  P2(Y2  +  z2f 

3/2 

(23) 


An  expression  for  the  sidewash  due  to  a  lifting  line  can  be  obtained  in 
a  similar  manner.  For  points  in  the  z  =  0  plane  not  directly  behind 
the  lifting  line,  equation  (23)  reduces  to 


0z(x,y,o) 


_i  fh2  (x  -  f)r(yi) 

2,r^h1  Y 2\/(x  -  f)2  -  p2Y2 


m 


For  points  directly  behind  the  lifting  line,  however,  the  equation  of  the 
lifting  line  must  be  known  in  order  to  evaluate  equation  (23)  in  the 
z  =  0  plane . 
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Unbent  lifting  lines.-  For  a  lifting  line  parallel  to  the  y-axis 
located  at  x  =  a,  equation  (23)  becomes 


,  f  rh2  ||x-  a)2-  P2Y2]y2-  z2[(x-a)2-  32(y2  +  2i2)]}(x-a)r 

jffa.rf.fr  Y2+22)]3/2  yi 


When  integrated  by  parts,  equation  (25)  becomes 


0z(x,y,z)  = 


XY(R2  -  P2z2)r(yi) 
2it(X2  -  (32z2)  (Y2  +  z2)R 


h2  XY(R2  -  02z2) 


•ar(yi) 


2,tJh1  (X2  -  p2z2)  (Y2  +  z2)R 

When  one  or  both  limits  of  the  first  term  of  equation  (26)  are  the  inter-  , 
section  of  the  Mach  cone  from  point  (x,y,z)  with  the  lifting  line,  these 
limits  are  neglected  because  only  the  finite  part  is  to  be  taken.  For 
lifting  lines  that  represent  airfoils,  T  is  zero  at  each  end  of  the  line; 
therefore,  in  this  case  equation  (26)  becomes 


rh2  XY(R2  -  P2z2)^il 

0Z(X>Y>Z)  =  -pi  /  — - - -  ^1 

Jhl  (X2  -  p2z2)(Y2  +  z2)R 


This  result  was  obtained  in  reference  d  by  using  vortex  theory 


dT(yj_)  zero,  equation  (26)  becomes 

dyx 


When 
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0z(x,y,z)  = 


xr 


{x2  -  p2  [(y  -  h2)2  +  2z2]}(y  -  h2) 


2rt(X  P  z  )  i  J(y  _  h2)2  +  z2j  ^2  _  p2jjy  -  h2)2  +  z2] 


{x2  -  P2 

(y  - 

-  hx)2  4  2z2 

(y  -  hi) 

[(y  -  hi)2  + 

■  z2 

\jx2  -  P2  [( 

y  -  hi)2  4  z2] 

(28) 


Equation  (28)  is  the  upwash  from  a  horseshoe  vortex;  this  result  was  first 
obtained  in  reference  11. 

The  velocity  potential  in  space  due  to  an  unbent  lifting  line  can  be 

dr 

expressed  in  terms  of  r  and  - .  Integrating  equation  (22)  by  parts 

ayi 

gives 


0(x,y,z)  =  -~y-1^  tan"1  — 
r  ’  2  Jt  YX 


dr 


h2 


2it 


>h2 


.  _l  zR  dT  , 
tan  YXd^dyi 


(29) 


For  a  horseshoe  vortex  -  is  zero  and 

dyi 


it hi)  =  r(h2) 

therefore,  the  potential  in  space  due  to  a  horseshoe  vortex  is 

„  _i  z\Ax  -  a)2  -  P2(y  -  h2)2  -  p2z2 

(y  -  h2)X 

_!  z\/(x  -  a)2  -  P2(y  -  hi)2  -  P2z2 

tan  x  — - - - - 

(y  -  hi)x 

The  components  of  the  velocity  can  he  found  from  equation  (30)  by 
differentiation. 


p 

0(x,y,z)  =  — 

2Jt 


(30) 


1 6 
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Lifting  line  of  constant  slope. -  The  potential  due  to  a  lifting  line  a 

of  constant  slope  can  he  expressed  as  (from  equation  (22)) 


where  the  equation  of  the  lifting  line  is 

y-]_  =  nu^  +  k 

When  integrated  hy  parts,  equation  (31)  becomes 
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The  potential  due  to  a  bent  lifting  line  as  shown  in  figure  5(a)  is 


r(o)t 


an 


-1 


mz 


V(X  +  e)2  -  p2(y2  +  z2) 

y(mx  +  k)  -  z2 


r(o)t 


an 


-i  mzV(x  "  g)2  ~ p2(y2  +  z2) 

y(mx  -  k)  -  z2 


no 


1_ 

2ir 


J-b/2 


dT 

dy- 


tan' 


a>\(x- 


■yi  +  k 


-  (32(Y2  +  z2) 


m 


dy- 


y(x  -  2111  ♦  !* 


m  /  m 


1 


>b/2 


z\/fx  -  — - ]  -  P2(Y2  +  z2) 

tan"1  — — - S—L - 


0 


ay  i 


dy 


y  x  - 


-  M 


'Oi.  /  v  m 


(33) 


18 


NACA  TN  2135 


where  m  is  the  slope  of  the  bent  lifting  line  when  y  is  positive 
(see  fig.  5)-  When  the  lifting  line  is  used  to  approximate  a  wing, 

IX-b/2)  and  T(b/2)  are  zero.  The  components  of  the  velocity  for  a  hent 
lifting  line  can  he  found  hy  differentiation  of  equation  (33)- 

When  d.n(yi)/dyi  is  zero  equation  (33)  becomes 


y(mx  +  k)  -  z^ 


Equation  (3*0  is  the  expression  for  the  potential  in  space  due  to  a 
vortex  of  constant  strength  as  shown  in  figure  5(h).  Since  the  velocity 
components  can  he  found  hy  differentiation,  the  upwash  is  given  hy 
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Approximations  to  the  Up wash  from  Lifting  Lines 

Although  the  integrals  that  arise  when  lifting  lines  are  used  are 
much  easier  to  evaluate  than  the  integrals  of  exact  lifting-surface 
expressions,  many  of  the  integrals  that  arise  in  connection  with  the 
use  of  lifting  lines  cannot  he  evaluated  in  closed  form  and  many 
involve  singularities .  It  is  therefore  desirable  to  obtain  approximate 
expressions  for  lifting  lines. 

Approximation  of  an  unbent  lifting  line.-  Equation  (28)  can  be  used 
to  approximate  an  unbent  lifting  line  by  a  series  of  horseshoe  vortices. 
The  upwash  from  a  series  of  horseshoe  vortices  that  approximate  the 
circulation  by  a  series  of  step  *  as  shown  in  figure  6(a)  is 


n 


i=0 


XYiQg  -  P2^2  -  2P2z2)[r(yi)  -  rfri-i)] 

(X2  -  p2z2)  (Yp2  +  z2)  \Jx 2  -  p^i2  -  P2z2 


(36) 


where  Yi  denotes  y  -  yj_  and  i  takes  on  all  integral  values  from  0 
to  n.  The  upwash  from  a  series  of  horseshoe  vortices  that  approximates 
the  circulation  by  a  series  of  steps  as  shown  in  figure  6(b)  is 


C  XY±(X2  -  pgYj2  -  2P2Z2)  [r(yi+i)  -  r(y±)] 
L—  (X2  -  p2z2)  (Yp2  +  z2)  \/x2  -  P2?!2  -  p2z2 


(37) 


The  average  of  equations  (36)  and  (37)  is 


r  XYj(x2  -  pgYj2  -  2pgz2)  [r(yi+1)  -  r(yp_iQ 
L—  (X2  -  P2z2) (Yp2  +  z2)  (/x2  -  p2(Yi2  +  z2) 


(38) 


Equation  (38)  is  the  upwash  from  a  series  of  horseshoe  vortices  which 
approximates  the  circulation  by  a  series  of  steps  as  shown  in  figure  6(c) 
In  general,  for  points  directly  behind  the  lifting  line,  equation  (38) 
should  give  best  agreement  with  the  exact  lifting  line  when  y  (the 
coordinate  for  the  field  point)  is  midway  between  consecutive  values  of 
y^.  An  expression  for  the  sidewash  may  easily  be  obtained  by  using  the 
sidewash  due  to  a  horseshoe  vortex  in  the  same  manner  that  the  upwash  due 
to  a  horseshoe  vortex  was  used  in  obtaining  equation  (38). 


-Approximation  of  bent  lifting  lines.  -  A  bent  lifting  line  can  be  approximated  by  a  series 
of  vortices  of  the  form  shown  in  figure  7  plus  terms  giving  the  effect  of  the  bend  in  the  center 
of  the  lifting  line.  Figure  8  shows  a  series  of  vortices  of  the  form  shown  in  figure  7  distrib¬ 
uted  along  a  bent  lifting  line.  The  upwash  from  such  a  system  of  vortices  is 
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-  3x/X  [(srast)  "  I)3Z  +  s(*  ”  ™)J  (sz  + 

(A  -  xtn)As2gEls  -  (giCgd  -  gx)  fsz  +  (A  -  xm)Aljm(o)j 
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where  ni  is  related  to  negative  values  of  y  and  n2  is  related  to 
positive  values  of  y.  The  first  term  on  the  right  side  represents  the 
vortices  for  negative  values  of  y;  the  second  represents  the  vortices 
for  positive  values  of  y.  The  remaining  terms  take  into  account  the 
effect  of  the  tend  in  the  center  of  the  lifting  line.  Equation  (39)  is 
set  up  for  a  distribution  of  vortices  as  shown  in  figure  8. 


Equation  (38)  can  be  used  to  approximate  a  bent  lifting  line  if  X 
is  replaced  by  where  X-^  denotes  x  -  xj..  When  this  is  done,  the 
result  is 


1 

'US 


n 


i=0 


XiYi(Xi2  -  P2*!2  -  2p2z2)  [r(yi+1)  -  rfy-^)] 
(Xi2  -  p2z2) (y^2  +  z2)  \^x i2  -  P2(Yi2  +  z2) 


(40) 


Equation  ( -U0 )  approximates  a  bent  lifting  line  by  adding  the  effects  of  a 
series  of  horseshoe  vortices  as  shown  in  figure  9(a)  when  the  circulation 
is  symmetrical  about  the  center  of  the  lifting  line,  and  when  the  points 
(xi>yi)  are  chosen  symmetrically  about  the  center  of  the  lifting  line. 

When  the  circulation  is  not  symmetrical  and/or  the  points  (xi,y^)  are 
not  chosen  symmetrically  about  the  center  of  the  lifting  line,  equation  (40) 
approximates  a  bent  lifting  line  by  a  series  of  vortices  as  shown  in 
figure  9(l>) •  The  downwash  from  an  infinite  line  vortex  of  constant 
strength  and  infinite  slope  is  zero;  therefore,  equation  (kO)  should  give 
a  reasonable  approximation  to  a  bent  lifting  line. 


Resume  and  Discussion  of  Theoretical  Development 

In  the  theoretical  development  exact  and  approximate  formulas  have 
been  developed  for  the  potential  and  the  upwash  due  to  thin  wings.  An 
expression  for  the  potential  in  space  due  to  a  thin  wing  was  derived  from 
an  expression  for  the  partial  derivative  of  the  potential  with  respect  to 
the  free-stream  direction.  The  upwash  was  found  by  differentiation  of  this 
expression.  Formulas  are  derived  for  the  upwash  in  the  wake  and  on  the 
plan  form.  The  results  of  the  exact  expression  for  thin  wings  were  used 
to  find  the  potential  and  the  upwash  due  to  lifting  lines,  and  the  exact 
expressions  for  lifting  lines  were  used  to  find  approximate  expressions 
for  lifting  lines. 

In  actual  calculations  better  approximations  to  lifting  surfaces 
may  be  obtained  by  using  several  lifting  lines  or  their  approximations 
distributed  along  the  chord.  Another  approach  would  be  to  distribute 
finite  vortices  over  the  plan  form.  Either  of  the  preceding  methods 
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may  be  used  to  find  a  good  approximation  to  the  downwash  from  a  lifting 
surface,  provided  that  the  apex  of  the  hyperbola  formed  by' the  inter¬ 
section  of  the  Mach  cone  with  the  plane  of  the  wing  does  not  lie  on  the 
plan  form.  For  cases  where  the  apex  is  on  the  plan  form  there  is  a 
finite  contribution  along  the  Mach  cone.  This  contribution  may  be 
evaluated  by  the  same  method  that  was  used  to  evaluate  the  double  inte¬ 
gral  in  equation  (15). 

It  should  be  remembered  that  this  paper  does  not  discuss  the 
effects  of  thickness;  however,  these  effects  may  be  evaluated  rather 
simply  by  use  of  source  distributions. 


APPLICATIONS 


The  results  of  the  theoretical  development  can  be  applied  to  a  number 
of  problems.  A  few  of  these  problems  will  be  considered  here. 


Potential  and  Upwash  at  Infinity 
When  x  approaches  infinity  equation  (6)  becomes 


0U,y,z)  =  1.  PP  ^X1  dxx  dyx 

2nJJ t  y2  +  z2 


T  Y2  +  z2 


Since 


£0(x,y,z)  =  f  A0X,  dx 
JL.E. 


integrating  the  preceding  expression  with  respect  to  gives 


«W>  -  f-  r  -2s-  ay 


hi  +  z2 


(M) 


Equation  ( hi)  is  the  same  equation  as  the  equation  for  the  potential  at 
infinity  behind  an  airfoil  in  subsonic  flow  and  is  a  well-known  result 
noted  in  references  1  and  2. 
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The  upwash  at  infinity  is  given  by 


(Y2  -  z2)A0ie 
(Y2  +  z2)2 


(42) 


Integrating  equation  (42)  by  parts  and  replacing  by  r  gives 


0z  = 


Yr 


2*(Y2 


+  Z‘ 


h2 


hl 


1_ 

2ir 


rh 2 


hi 


dT 

3yi 


y2  +  z2 


dyi 


(43) 


For  wings  r(hi)  and  r(h2)  are  zero,  the  first  term  of  equation  (43) 
becomes  zero,  and  the  remaining  term  gives  the  upwash  at  infinity.  The 
same  result  is  given  in  reference  6. 


¥ 


The  Upwash  for  a  Flat  Unswept  Wing  of  Infinite  Span 

*2n\T 

The  A0X  on  the  plan  form  of  an  infinite  unswept  wing  is  — — . 

Since 

dA0x]_ 

- - ±  =  0 

oyi 

equation  (17)  reduces  to 

0z(x>y>°)  =  "aV 

For  points  in  the  wake,  equation  (14)  gives 

0z(x>y>°)  =  0 

The  preceding  results  are  the  well-known  equations  for  the  upwash  in  the 
2=0  plane  for  the  two-dimensional,  flat,  unswept  wing. 

t 

Camber  and  Twist  Necessary  to  Give  Specified 
Lift  Distributions 

Equation  (17)  may  be  used  to  find  the  camber  and  twist  necessary  to 
give  a  specified  lift  distribution.  The  local  angle  of  attack  is  given 
by 

-  A 

V 


cr 


(wo 


NACA  TN  2135 


25 


\ 


and  the  pressure  coefficient  is  given  by 


P 


2A0xi 

V 


(^5) 


The  local  angle  of  attack  for  a  given  lift  distribution,  derived  from 
equations  (17) ,  (44),  and  (45),  is 


m 


As  an  example,  the  local  angle  of  attack  of  a  rectangular  wing  with 
constant  load  is  now  considered. 

When  the  pressure  distribution  is  a  constant,  equation  (46)  becomes 


For 

equation 


the  rectangular  wing  in  the  region  unaffected  by  the  tip, 
(47)  becomes 


■  -  -! 


(48) 


In  the  region  affected  by  the  tip,  equation  (47)  becomes 

flO  - - : - - 

\j(x  -  xi)2  -  P2^2 


cr  = 


k  kity 


x+py 


(x  -  xx) 


dx- 


(^9) 


where  the  coordinates  are  now  located  at  the  intersection  of  the  leading 
edge  and  the  tip.  (See  fig.  10.)  After  the  integration  is  performed 
equation  (49)  becomes 


Equation  (50)  agrees  with  the  results  in  reference  10. 
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Upwash  Close  to  the  Trailing  Edge 

The  value  of  the  vertical  perturbation  velocity  immediately  behind 
the  trailing  edge  can  be  obtained  from  equation  (17). 

In  this  application  the  concept  of  cancellation  of  pressure  is 
utilized.  It  is  assumed  that  the  wing  extends  past  the  trailing  edge. 
Since  the  pressure  must  be  zero  through  the  wake,  this  pressure  behind 
the  trailing  edge  must  be  canceled.  The  boundary  conditions  in  the  wake 
are  then  satisfied  by  the  wing  pressure  extended  past  the  trailing  edge 
and  the  cancellation  pressure.  The  upwash  is  therefore  made  up  of  the 
effect  of  the  pressure  on  the  wing,  its  extension  behind  the  wing,  and 
the  cancellation  of  the  pressure  behind  the  trailing  edge. 

For  trailing  edges  which  are  perpendicular  to  the  free-stream  direc¬ 
tion  the  application  of  equation  (17)  (using  the  cancellation-of -pressure 
concept)  to  the  upwash  immediately  behind  the  trailing  edge  leads  directly 
to  the  relation 


(0z)ije  -  -ctteV  +  2^x^TE 


(51) 


./here  is  the  local  angle  of  attack  of  the  wing  at  the  trailing  edge. 

Che  result  given  by  equation  (51)  was  obtained  in  references  1  and  2  from 
1  physical  consideration.  When  the  trailing  edge  has  a  slope  m,  equation 
(17)  yields 


(0z)ipE 


-ff^V  + 
TE 


(a0x)te  \/p21112  -  1 
“  2m 


(52) 


Squation  (52)  is  also  obtained  in  reference  1. 


Upwash  in  the  Wake  behind  a  Rectangular  Wing 
Pitching  about  Its  Leading  Edge 

As  an  example  to  illustrate  downwash  calculations  the  upwash  in  the 
plane  of  the  wing  from  a  rectangular -pitching- wing  is  now  calculated. 
(See  fig.  10  for  axes  used  in  analysis.)  For  stationary  axes  a  pitching 
wing  moves  in  the  arc  of  a ’circle;  therefore,  for  axes  fixed  in  the  wing 
a  pitching  wing  in  linearized  flow  can  be  replaced  by  a  wing  which  has  a 
local  angle  of  attack  that  varies  linearly  in  the  free-stream  direction. 
It  is  assumed  that  the  wake  remains  in  the  z  =  0  plane  and  that  the 
rolling-up  of  the  trailing  vortices  can  be  neglected*. 
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The  wake  behind  a  rectangular  wing  may  be  divided  into  two  regions. 
One  region  is  not  affected  by  the  tip  or  tips  and  is  therefore  a  region 
which  has  a  two-dimensional  flow.  The  other  region  is  the  region  affected 
fcy  tip  or  tips;  in  this  region  the  flow  is  three-dimensional. 

There  is  no  upwash  in  the  wake  behind  a  two-dimensional  unswept, 
pitching  wing;  therefore,  only  the  region  affected  by  the  tips  of  the 
wing  has  upwash.  Since  the  analysis  is  based  on  linearized  flow,  the 
two  tip  effects  add  directly  and  therefore  only  one  tip  need  be  treated. 

The  upwash  in  the  z  =  0  plane  due  to  one  tip  of  a  rectangular 
pitching  wing  is  shown  in  figures  10  and  11.  The  velocity  potential  and 
the  pressure  on  a  rectangular  wing  pitching  about  the  leading  edge  were 
obtained  by  transforming  the  expressions  given  in  reference  12  for  the 
velocity  potential  and  the  pressure  on  a  rectangular  wing  pitching  about 
the  half-chord  line. 

The  upwash  at  the  trailing  edge  was  found  by  using  equation  (51)* 

The  upwash  close  to  the  trailing  edge  was  found  by  use  of  equation  (17). 

The  evaluations  of  the  integrals  were  performed  numerically.  Equation  (17) 
was  derived  to  find  the  downwash  on  the  plan  form,  but  it  may  be  used 
in  the  wake  by  applying  the  cancellation-of -pressure  concept. 

A  lifting  line  was  used  for  calculating  the  upwash  2  or  more  chords 
behind  the  trailing  edge  inasmuch  as  this  method  was  found  to  be  accurate 
to  two  decimal  places  in  this  region. 

The  upwash  at  infinity  was  found  by  using  equation  (43).  At  infinity 
the  upwash  in  the  z  =  0  plane  for  0  =  C  =  1  is  given  by  the  following 
expressions: 

For  y  >  0 


(53a) 


for  -1  <  y  <  0 


(53b) 


for  y  <  -1 
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Approximation  of  the  Upwash  from  a 
Rectangular  Pitching  Wing 

The  calculation  of  the  downwash  from  the  rectangular  pitching  wing 
together  with  results  of  reference  6  indicates  that,  for  regions  some 
distance  behind  the  trailing  edge,  lifting  lines  give  very  good  approxi¬ 
mations  to  the  actual  flow  fields.  In  many  cases  it  is  simpler  to  use 
the  approximation  given  by  equation  (38)  "than  to  use  the  exact  lifting¬ 
line  expression.  The  accuracy  of  equation  (38)  was  investigated  by  cal¬ 
culating  the  upwash  from  an  unbent  lifting  line  and  comparing  the  results 
with  the  results  obtained  using  equation  .(38).  The  unbent  lifting  line 
used  was  one  that  approximates  the  flow  from  the  pitching  rectangular 
wing. 

In  figure  12  the  values  of  the  upwash  from  the  exact  and  the  approxi¬ 
mate  lifting  lines  are  plotted.  Ten  equally  spaced  points  across  the  tip 
region  were  used  in  the  approximation.  It  can  be  seen  that  for  this  case 
in  most  regions  the  approximation  yields  results  which  are  as  reliable  as 
the  exact  values.  If  greater  accuracy  is  required,  it  should  be  remem¬ 
bered  that  the  accuracy  of  the  approximation  increases  as  the  number  of 
points  is  increased. 

Approximation  to  the  Upwash  from  a 
Bent  Lifting  Line 

The  approximation  to  the  upwash  from  a  bent  lifting  line  given  by 
equation'  (39)  was  used  to  find  the  downwash  along  the  line  z  =  0,  y  =  0 
behind  a  triangular  wing  as  shown  in  figure  13.  Figure  13  shows  the 
values  given  by  equation  (39)  for  10  points  across  the  semispan,  the  exact 
lifting  line,  and  the  exact  linearized  solution.  The  values  for  the  exact 
lifting  line  for  this  case  were  taken  from  reference  d  and  the  values  of 
the  exact  linearized  solution  were  taken  from  reference  3. 

Figure  13  shows  good  agreement  between  the  approximate  and  exact 
bent  lifting  lines.  In  general,  the  accuracy  of  the  approximation  will 
increase  as  the  number  of  points  used  in  the  approximation  is  increased; 
therefore,  it  is  expected  that  the  difference  between  the  exact  and 
approximate  expressions  given  in  figure  13  could  be  decreased  by  increas¬ 
ing  the  number  of  points  across  the  semispan. 
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APPENDIX 

DIFFERENTIATION  UNDER  THE  INTEGRAL  SIGN 
OF  THE  EXPRESSION  FOR  0 

The  fact  that  equation  (6)  may  be  differentiated  with  respect  to  z 
under  the  integral  sign  without  consideration  of  the  variable  limits  can 
be  proved  as  follows: 

If  it  is  assumed  that  the  right  side  of  equation  (6)  can  be 
differentiated  by  differentiating  only  the  integrand,  then 


(Al) 


If  the  preceding  equation  is  correct,  the  potential  may  be  expressed  by 
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The  result  of  performing  the  first  integration  in  equation  (A3)  is 


Equation  (A4)  is  the  same  as  equation  (6);  therefore,  it  has  been  shown 
that  equation  (6)  may  he  differentiated  with  respect  to  z  under  the 
integral  sign  without  regard  to  the  variable  limits.  Similar  proof  for 
the  other  velocity  components  can  also  be  obtained. 
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Mach  cone 


Intersection  of  Mach  cone 
with  plan  form 


plan  form 
plane 


y,=y 


(a)  Arbitrary  point  above  plan  form.  61^.105  1 

Figure  1.-  The  intersection  of  forward  Mach  cone  from  point  (x,y, z)  with 

arbitrary  plan  form. 
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(a)  Bent  lifting  line. 


(b)  Vortex  of  constant  strength. 
Figure  5.-  Lifting  line  of  constant  slope. 
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(a)  Approximation  that  leads  to  equation  (36). 


(t>)  Approximation  that  leads  to  equation  (37). 


.  (c)  Approximation  that  leads  to  equation  (38). 

Figure  6.-  Approximations  of  the  span  loading  by  series  of  rectangles. 
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(b)  Unsymmetrical  loading  and/or  unsymmetrically  distributed 

points . 


Figure  9.-  Bent  lifting  lines  approximated  by  a  series  of  horseshoe 

vortices . 
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Figure  12.-  Downwash  from  a  lifting  line  and  its  approximation  by  a  series 
of  horseshoe  vortices  at  ten  equally  spaced  points.  (Lifting  line 
approximates  a  rectangular  pitching  wing.) 


Figure  13.-  Downwash  determined  by  exact  and  approximate  lifting  lines 
alone  the  line  z  -  0,  y  =  0  behind  a  triangular  wing  with  aspect 
ratio  of  3.2.  M  =  \/2. 
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